The Fermi problem with artificial atoms in circuit QED 
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We propose a feasible experimental test of a 1-D version of the Fermi problem using supercon- 
ducting qubits. We give an explicit non-perturbative proof of strict causality in this model, showing 
that the probability of excitation of a two-level artificial atom with a dipolar coupling to a quantum 
field is completely independent of the other qubit until signals from it may arrive. We explain why 
this is in perfect agreement with the existence of nonlocal correlations and previous results which 
were used to claim apparent causality problems for Fermi's two-atom system. 
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Information cannot travel faster than light. But in 
quantum theory, correlations may be established between 
spacelike separated events. These facts are not contradic- 
tory, since correlations need to be assisted with classical 
communication in order to transmit information. 

The two physical phenomena above arise in a natural 
fashion in the following situation, which is the so-called 
Fermi problem [1 j, originally proposed by Fermi to check 
causality at a microscopic level. At t = a two-level neu- 
tral atom A is in its excited state and a two-level neutral 
atom B in its ground state, with no photons present. If A 
and B are separated by a distance r and v is the speed of 
light, can A excite B at times t < r/v? Fermi 's answer 
was negative but his argument had a mathematical flaw. 
When a proper analysis is carried on, fundamental quan- 
tum theory questions arise due to the interplay between 
causal signaling and quantum non-local phenomena. 

These issues led to a controversy [2H5] on the causal 
behavior of the excitation probability of qubit B, whose 
conclusions were never put to experimental test. A no- 
torious claim on causality problems in Fermi's two-atom 
system was given in [2]. The reply of [3 was in the ab- 
stract language of algebraic field theory and the proof of 
strict causality of [5 is perturbative, although given to 
all orders in perturbation theory. The Fermi problem is 
usually regarded just as a gedanken experiment, and re- 
mains untested, essentially because interactions between 
real atoms cannot be switched on and off. 

With this work we give a complete description of the 
problem in a physical framework in which predictions can 
be verified. This framework will be circuit QED which 
can be regarded as a 1-D version of Quantum Electrody- 
namics (QED) with two-level (artificial) atoms, a testbed 
which makes it possible to control the interaction and 
tune the physical parameters. We complete previous de- 
scriptions made of the problem and explain how there are 
no real causality issues for Fermi's two-atom system. We 
give an explicit non-perturbative proof of strict causality 
in these setups, showing that the probability of excita- 
tion of qubit B is completely independent of qubit A 
for times t < r/v and for arbitrary initial states. As 
a matter of fact, this comes as a manifestation of the 



nonsignaling character of the quantum theory [6]. We 
also show how this is compatible with the existence of 
nonlocal correlations at times 0<t<r/v,a fact pointed 
out in various theoretical proposals to entangle qubits at 
arbitrarily short times [TUTU]. More precisely, we give a 
non-perturbative proof of the fact that the probability of 
B being excited and A in the ground state is finite and in- 
dependent at any time, even for t < r/v. We provide also 
a physical and intuitive explanation of why the conclu- 
sions in [2], even if mathematically sound, do not apply 
to the causality problem. At the end we discuss the time 
dependence predicted in our model for the various exci- 
tation probabilities and suggest a feasible experimental 
test of causality using superconducting circuits. 

In what follows we focus on a practical setup of circuit- 
QED, with two qubits, A and £?, interacting via a quan- 
tum field. The qubits have two stationary states \e) and 
\g) separated by an energy HQ and interact with a one- 
dimensional field, V(x), which propagates along a one- 
dimensional microwave guide that connects them 

/oo 
dk ^Nuo k e ikx a k + H.c. (1) 
-oo 

This field has a continuum of Fock operators [afc,aj,,] = 
5(k — k'), and a linear spectrum, uj k = v\k\, where v is the 
propagation velocity of the field. The normalization and 
the speed of photons, v = (c/) _1//2 , depend on the micro- 
scopic details such as the capacitance and inductance per 
unit length, c and I. We will assume qubits that are much 
smaller than the relevant wavelengths, A = v/Q, and are 
well separated. Under these conditions the Hamiltonian, 
H = Ho + Hi, splits into a free part for the qubits and 
the field 

1 f°° 
H = -h£t{(7 z A + cf z b ) + / dk fujj k a\a k (2) 

^ J -oo 

and a point-like interaction between them 

Hi= Yl dja x jV{xj) (3) 

J=A,B 

Here xa and xb are the fixed positions of the atoms, and 
dj g x j is equivalent to the dipole moment in the case of 
atoms interacting with the electromagnetic field. 
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The original formulation of the Fermi problem begins 
with an initial state 



\in) 



(4) 



in which only qubit A has been excited, while B and the 
field remain in their ground and vacuum states, respec- 
tively. The total probability of excitation of qubit J is 
the expectation value of the projector onto the excited 
state Vj = \ej) (ej \ . In the Heisenberg picture 



P eJ = (in\V e j(t)\in), Je{A,B}. 



(5) 



We will prove that for vt < r the probability P e B is com- 
pletely independent of the state of qubit A for all initial 
states. In the Heisenberg picture this amounts to show- 
ing that there appears no observable of A in the projector 
V B (t) for vt < r. Our proof begins by solving formally 
the Heisenberg equations for Vj 



V e j{t) - V e j(0) = 



dj 



I dt'a V j(t')V(xj,t'). 
Jo 



(6) 



Integrating the Heisenberg equations of the operators a k 
and a\ and inserting them in Eq. the total field eval- 
uated at x in Heisenberg picture is decomposed 



V(x,t) = V (x,t) + V A {x,t) + V B (x,t) 
into the homogenous part of the field 

V (x,t) = 



H.c. 



(7) 



(8) 



and the back-action of A and B onto the field 
-idjN 



Vj(x,t) = 



h 



JO J-i 



(9) 

Uke ik(x-xj)-iu> k (t-t') dkdt > + H c 



Eqs. ^ translates into a similar decomposition for the 
probabilty V% with three terms 



n(t) = n (t)+n B (t)+nA(t) 



(10) 



which are proportional to Vo,Vb and Va, respectively. 
The only explicit dependence on A may come from V BA 
through Va Manipulating Eq. ([9| gives 



V A {x B ,t) 



-2nd A N d 
h dr 



(ii) 



where the Heaviside function shows that strictly 
T^ BA (x B ^t) = for vt < r, and no such dependence is 
possible. We still have to analyze a possible implicit de- 
pendence on A through V BB , whose expression is 



r e BB (t) 



id 2 B N 



dt' 



dt"o- y B {t')o-Ut") 
dkuj k e- iu)h(t, - tff) +H.c. 



(12) 
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Figure 1: The terms of order dj,d 2 j and d 3 j contributing to 
the amplitude for exciting qubit B. 



The only implicit dependence could come through the 
evolution of a B y (t), but again this is not the case. Since 
[a B ,Hi] = 0, the evolution of a B does not involve the 
field in any way, and for cr B (t) we have that & 1 L(t) = 
^V(x Bl t)a B (t) so using again Eq. (M) and 

BB 



Qa B (t)/2 

Eq. (11) we see that the A-dependent part of f BB is 
for vt < r. Thus T B may be finite but is completely 
independent of qubit A for vt < r, as we wanted to show. 

So far, we have demonstrated that although P B (t) is 
non-zero for vt < r, the only non-zero contribution is 
V B0 , which is not sensitive to the qubit A and thus can- 
not be used to transmit information between the qubit s. 
Now we will show that this result is compatible with the 
existence of correlations for vt < r. For instance, we con- 
sider the probability of finding qubit B excited and qubit 
A on the ground state P e B,gA, which is: 



eB,gA 



{in\V%{t)V 9 A (t)\in 



(13) 



where V g A 



1 



V%. Using Eqs. Q and (To]), we 



find a term in this probability which is proportional to 
V e BB V 9 AA and thus to V B {x B ,t) V A (x A ,i). From Eq. ^ 
we obtain: Vj(xj,t) ex ^{o~j(t)6(t)}. Therefore, we con- 
clude that in ( [13] ) there is an unavoidable dependence on 
A at any t > 0, but this is not a causality violation be- 
cause correlations alone cannot transmit information. 

At this point it remains a single question: How can the 
A-dependent part of P eB be zero while the one of P e B, g A 
is nonzero for vt < r? To better understand it we need 
less formal results that rely on perturbative expansions, 
but we would like to remark here that the conclusions 
above are valid to all orders in perturbation theory. To 
obtain the total probability of excitation of qubit B P eB 
to a given order in perturbation theory, one has to expand 
to a certain order the operators appearing in Eqs. (|6|, 
([7]), (JoJ) . The different terms in the expansion can be 
related to the probabilities of the different physical pro- 
cesses involved. Fig.[T]shows the diagrams of the different 
amplitudes contributing to P eB up to the fourth order in 
dj. The lowest order amplitude contributing to a final 
excited B qubit is of order dj, which means that terms 
up to order dj have to be considered. The only terms 
leading to this final state will be = v B , = 



3 



x + ua vb, = cl'vb + UaVaVb + V B UbVb + <LW^ 

where uj (vj) represent the amplitude for single pho- 
ton emission at qubit J when the qubit is initially in 
the ground (excited) state, x is the amplitude for pho- 
ton exchange, aj are the radiative corrections of qubit J, 
and finally 5A4^ is the amplitude for photon exchange 
accompanied by a single photon emission at qubit A. No- 
tice that some of these processes are only possible beyond 
the rotating wave approximation, which breaks down for 
strongly coupled circuit- QED setups [18] as the ones con- 
sidered later. Keeping only terms up to fourth order, we 
have for the probability to get B excited at a time t 



PeB(t) 



\M^\ 2 + \M^\ 2 +2Re{M^*MW} 
2Re{M {ir M {3) } + 0(d 5 ) (14) 



The final states in are orthogonal to those in Ai^ 

and to the three photon terms in A^ 3 ). Hence, their 
interference vanishes. Besides, we are only interested in 
the A-dependent part of the probability, so we can re- 



move the r-independent terms left in (14), marking the 
remaining contributions with a superscript ( r ) 



9 (r) 



-2Re{M iir SM {3) } + (D(d 5 ). (15) 



The first term actually gives P e B, g A up to the fourth or- 
der, it is positive and A-dependent at all times, as shown 
in Fig. [2^l. The second term is not a projector onto any 
physical state, but an interference term which has the 
effect of canceling out exactly the first term for vt < r 
but not for vt > r (cf. Fig. |2|d). In a nutshell, interfer- 
ence seems to be the physical mechanism that operates 
at all orders in perturbation theory to give the causal be- 
havior of the total probability of excitation that we had 
previously shown. 

These perturbative results cast new light on the con- 
troversy on the Fermi problem and help us understand 
why our results do not contradict those of Hegerfeldt [2] . 
Hegerfeldt proved mathematically that the expectation 
value of an operator consisting of a sum of projectors 
cannot be zero for all the times vt < r, unless it is zero at 
any time. Indeed, the expectation value of V%{t) cannot 
be zero for all vt < r, for it always contains the contri- 
bution Vbo from Eq. (JsJ) - However, as we showed non- 
perturbatively, the actual relevant question for causality 
is whether the expectation value of V% A (t) vanishes for 
vt < r or not, since only this part of the probability is sen- 
sitive to qubit A and could be used to transmit informa- 
tion. Besides, according to our above perturbative results 
to fourth order, the r-dependent part of the probability, 
that is the expectation value of V% A (t), is not a mere sum 
of projectors, but also contains interfering terms. Thus, 
Hegerfeldt's result does not apply and V% A (t) can ^ e zero 
for vt < r as is actually the case. Both results are in ac- 
cord with a general fact of Relativistic Quantum Field 
Theory: two global states can not be distinguished lo- 
cally with the aid of a local projector annihilating one of 




Figure 2: (a) P^ A and (b) P^ versus vt/r for Qr/v = | 
(blue, crosses), tt (red, squares), and 2tt (green, circles) with 
Ka,b = 0.0225. For vt < r the qubits are spacelike separated, 
but there are correlations between them and figure (b) shows 
the expected causal behavior, (c) P^ (blue, circles) and 
|7W (1) | 2 (green, squares) vs. vt/r for K A = 0.20, K B = 0.04 
and a separation of one wavelength r = 2tvv/Qa,b- With this 
data and £1/2tt c± 1GHz we have At ~ Ins. (Color online) 



(r) 



the states, since the local observable algebras are Type 
III von Neumann algebras (See [3[ H] for a discussion) . 
We will now suggest an experiment to check the causal 
behavior of P e s- For this we need to control the inter- 
action time at will to access the regions at both sides 
of t = r/v. This, which is highly unrealistic with real 
atoms, becomes feasible in circuit-QED. While the ideas 
are valid for both inductive and capacitive couplings, we 
will focus on using a pair of three-j unction flux qubits 
[TT| H2] . Each of the qubits is governed by the Hamilto- 
nian H 0J = \ejG z j + \Ajo x j. The energy ej = 2I p 5$ xJ , is 
approximately linear in the external magnetic flux, 5& x j, 
measured from the degeneracy point, and we assume that 
the gap Aj is fixed. The result is a qubit energy differ- 
ence nj(5$ xJ ) = v /(2/p^,j) 2 + A 2 J . 

The coupling between the qubit and the microwave 
photons is ruled by the dimensionless ratio 



K.j 



-^-2(5/0/) 



(16) 



Here g = djVNft/H is the coupling strength between 
a qubit and the cavity that would be obtained by cut- 
ting the transmission line to be perfectly resonant with 
the qubit transition. These numbers enter the qubit ex- 



citation probability computed before (15) through the 
product P^{t) oc K A K B . Since Kj oc l/f2j, we may 
use the external fluxes to move from a weakly coupled 
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regime with no qubit excitations, Qj <C g, to the maxi- 
mum coupling strength, Qj ~ Aj ((S^j = 0). 

Let us first discuss how to prepare the initial state 
Q of the Fermi problem. We assume that the system 
starts in a ground state of the form \gAgB®)- This is 
achieved cooling with a large negative value of 5$> x j 
on both qubits, which ensures a small value of g/Vtj 
and Kj. We estimate that couplings g/flj < 0.15 and 
Qj ~ 1.5GHz lower the probability of finding photons 
in the initial state below 5 x 10 -3 , both for vacuum and 
thermal excitations. Both magnetic fluxes are then raised 
up linearly in time, 5& x j = ajt, to prepare the qubits. 
Using a Landau-Zener analysis [13 of the process we con- 
clude that an adiabatic ramp q;^ <C nA B /4hl p of qubit 
B followed by a diabatic ramp p~4j [15] a a ^> A\/h2I p 
of qubit A, leads to the desired state \eAgB 0) with a fi- 
delity that can be close to 1, depending only of a^, as 
derived from the Landau-Zener formula. Note that the 
minimum gap Ab has to be large enough to ensure that 
the qubit-line coupling of B remains weak and the qubit 
does not "dress" the field with photons. 

Once we have the initial state, both magnetic fluxes 
must take a constant value during the desired interac- 
tion time. After that, measurements of the probability 
of excitation of qubit B can be performed with a pulsed 
DC-SQUID scheme [EE [17]. The timescale of the "jump" 
of the probability around t = r/v for qubit frequencies 
in the range of GHz and a separation of one wavelength 
r = 27rv/n is At ~ Ins [Fig. ^} . Although the total 
measurement of the SQUID may take a few /is, the cru- 
cial part is the activation pulse (~ 15ns) in which the 
SQUID approaches its critical current and may switch 
depending on the qubit state. During this activation pe- 
riod the SQUID and the qubit are very strongly coupled 
(g ~ GHz), [16] and the qubit is effectively frozen. The 
time resolution of the measurement is thus determined by 
the ramp time of the activation pulse, which may be be- 
low nanoseconds. Among the sources of noise that are ex- 
pected, the short duration of the experiment, well below 
Ti and T2 of usual qubits, makes the ambient noise and 
decoherence pretty much irrelevant. Thermal excitations 
of the qubits and the line may be strongly suppressed by 
using larger frequencies (> 1.5GHz). The most challeng- 
ing aspect is the low accuracy of SQUID measurements, 
which are stochastic, have moderate visibilities [16] and 
will demand a large and careful statistics. 

On the technical side, it is important to choose care- 
fully the coupling regimes. If we wish to compare with 
perturbation theory, we need Kj <C 1. However, at the 
same time the product KaKb must take sizable values 
for P e B oc KaKb to be large. And we need to dis- 
criminate the causal signal from the r-independent back- 
ground of the probability of excitation, whose main con- 
tribution is lA^ 1 ^ 2 oc Kb- Thus, a good strategy would 
be to work with Ka > Kb- In Fig. ^ we show that it 
is possible to achieve a regime in which the perturbative 



approximations are still valid and the r-dependent part 
of P c b is comparable to |A^^^| 2 in the spacetime region 
of interest vt ~r. 

To conclude, we have considered a system of two super- 
conducting qubits coupled to a transmission line, which 
can be suitably described in the framework of 1-D QED 
with two- level (artificial) atoms. Starting from an initial 
state with qubit A excited, qubit B in the ground state 
and no photons, we have illustrated the causal character 
of the model showing that the probability of excitation 
of qubit B is completely independent of qubit A when 
vt < r. We have also shown that this is in agreement 
with the existence of nonlocal correlations and we have 
used perturbative computations to see the physical mech- 
anism underlying the causal behavior. Finally, we have 
suggested an experiment feasible with current technology 
that would solve the controversy on the Fermi problem. 
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